Dynamics and Thermodynamics of the Low- Temperature 
Strongly Interacting Bose Gas 
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We measure the zero-temperature equation of state of a homogeneous Bose gas of ^Li atoms by 
analyzing the m-situ density distributions of trapped samples. For increasing repulsive interactions 
our data shows a clear departure from mean-field theory and provides a quantitative test of the many- 
body corrections first predicted in 1957 by Lee, Huang and Yang. We further probe the dynamic 
response of the Bose gas to a varying interaction strength and compare it to simple theoretical 
models. We deduce a lower bound for the value of the universal constant ^ > 0.44(8) that would 
characterize the universal Bose gas at the unitary limit. 
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From sand piles to neuronal networks, electrons in 
metals and quantum liquids, one of the greatest chal- 
lenges in modern physics is to understand the behav- 
ior of strongly interacting systems. Recent advances in 
quantum gas experiments have made it possible to test 
with unprecedented accuracy complex many-body theo- 
ries beyond a simple mean-field description. Strong cor- 
relations in quantum gases can be achieved for example 
by loading the atoms into an optical lattice, where the 
enhanced role of the interaction with respect to the ki- 
netic energy can drive the gas into a strongly correlated 
Mott insulating state [1]. Another possibility is to re- 
duce the dimensionality to increase quantum fluctuations 
in the system [2|. In the fundamental case of a homo- 
geneous three-dimensional Bose g famous example 
is superfluid "^He, where the short-range interaction po- 
tential plays a crucial role. For dilute ultracold atomic 
vapors, the properties of the strongly interacting Bose 
gas remain largely to be explored. Magnetic Feshbach 
resonances make it possible to tune the atom-atom in- 
teractions, characterized by the s-wave scattering length 
a However, first studies of the strongly interact- 

ing regime were plagued by severe inelastic atom loss 
for Bose systems [5|. In the last two years, experiments 
at JILA and Rice have revived interest in these systems 
0,0]. Here we report on a quantitative measurement of 
the thermodynamic Equation of State (EoS) of a strongly 
interacting atomic Bose gas in the low-temperature limit. 
In a first part, we restrict ourselves to a regime with neg- 
ligible atom loss {a/ao < 2150, oq being the Bohr radius), 
and our EoS reveals the first beyond mean-field correc- 
tion due to quantum fluctuations. We perform Quantum 



Monte Carlo (QMC) simulations to support our experi- 
mental flndings. In a second part we discuss the dynam- 
ics of a Bose gas brought to the unitary regime. 
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FIG. 1: (Color online) Doubly-integrated density profile of a 
trapped Bose gas at a scattering length a/ao ~ 2150. The 
average over 5 experimental images is shown in black points. 
The QMC predictions for 3.9 10* atoms are plotted in solid 
line for T/Tc = 0.75 in red, 0.5 in orange, 0.25 in green and 
0.125 in purple. Inset: deviation per degree of freedom 
of a single experimental density profile with QMC results at 
different temperatures. The excellent agreement between ex- 
perimental profiles and QMC validates the zero-temperature 
assumption for the EoS measurement. 

A seminal result for the thermodynamics of the dilute 
Bose gas is the expansion of the ground state energy (per 
volume V) in terms of the canonical gas parameter na^, 
first obtained in the late '50s Q: 



(1) 



where n is the density of the gas, and g — Anh a/m 
the coupling constant for particles with mass m. The 
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FIG. 2: (Color online) Equation of State of the homogeneous 
Bose gas expressed as the normalized pressure h as a function 
of the gas parameter v. The gas samples for the data shown 
in the main panel (resp. upper inset) have been prepared at 
scattering lengths of a/ao ~ 1450 and 2150 (resp. a/ao — 
700). The red sohd line corresponds to the LHY prediction, 
and the red dashed line, to the mean-field EoS h{u) — 2iny^. 
The solid black line represents the QMC EoS at T/Tc = 0.25 
and is nearly indistinguishable from the LHY EoS (see lower 
inset). The green area delimits the uncertainty of 5% on the 
value of a. 



first term in EqUJis the mean- field energy, while the Lee- 
Huang- Yang (LHY) correction, proportional to \/ na^, is 
due to quantum fluctuations Q. Up to this order, the ex- 
pansion is universal^ in the sense that it solely depends 
on the gas parameter na^ and not on microscopic de- 
tails of the interaction potential 0-ll[. The term be- 
yond LHY involves the first non-universal correction 12 1. 
Despite its fundamental importance, this expansion has 
never been checked experimentally for an atomic Bose 
gas. 

Our experimental setup was described in [l3| . Starting 
from a ^Li cloud in a magneto-optical trap, we optically 
pump the atoms into the \F — 2,mp = 2) hyperfine state 
and transfer them into a magnetic loffe trap. After evap- 
orative cooling to a temperature of ~ 4 /iK, the atoms 
are loaded into a hybrid magnetic/optical trap and then 
transferred to the \F — 1, mp — 1) state. The radial op- 
tical confinement of the trap is provided by a single laser 
beam of waist 35 fXTn operating at a wavelength 1073 nni, 
while the weak axial confinement is enhanced by an ad- 
ditional magnetic curvature. We apply a homogeneous 
magnetic field to tune the interaction strength by means 
of a wide Feshbach resonance around 740 G. The final 
stage of evaporation in the optical trap is carried out at 
a bias field of 717 G, where the scattering length has a 
value of about 200 ao, and results in a Bose- Einstein con- 
densate of ~ 6 X 10^ atoms with no discernible thermal 
part. In the final configuration the trapping frequen- 
cies are given by = 27r x 345(20) Hz in the radial and 
ljJz = 27r X 18.5(1) Hz in the axial direction. The magnetic 
bias field is then adiabatically ramped to the vicinity of 



the Feshbach resonance and the density distribution is 
recorded using in-situ absorption imaging. 

For the measurement of the EoS, we follow the method 
of 14- T^l- Accordingly, the local pressure P{z) along the 
symmetry axis of a harmonically trapped gas is related 
to the doubly-integrated in-situ density profile n{z) — 
J dx dy n{x, y, z): 



(2) 



This formula relies on the local-density approximation 
in which the local chemical potential is defined as Hz = 
^0 — V{0,0, z), where /xq is the global chemical potential 
of the gas, and V^(r) is the harmonic trapping potential. 

As the EoS critically depends on the scattering length, 
a precise knowledge of the latter close to the Feshbach 
resonance of the |1,1) state is essential. In view of the 
discrepancy between two recent works 0, [l9| . we have 
independently calibrated the scattering length a{B) as a 
function of magnetic field B by radio-frequency molecule 
association spectroscopy [l^ , as described in the supple- 
mentary material. To measure the pressure at different 
interaction strengths we have selected images with atom 
numbers in the range of 3 - 4 x 10^ in order to avoid high 
optical densities during absorption imaging while keep- 
ing a good signal-to-noise ratio. A total of 50 images are 
used, spanning values of a/ao from 700 to 2150. We cali- 
brate the relation between the integrated optical density 
and the pressure of the gas at weak interaction, well de- 
scribed by mean-field theory (inset of Figl2]). The density 
profiles then generate the EoS ([2]). The global chemical 
potential /^o remains to be determined. For this work, 
we infer /io self-consistently in a model-independent way 
from the density profiles (see the supplementary mate- 
rial) . 

In the dilute limit na^ <C 1, where the EoS is universal, 
the grand-canonical EoS of the homogeneous Bose gas at 
zero temperature can dimensionally be written as 



Pi^i,a) 



■h{iy=^a^ 
9 



(3) 



where v is the (grand-canonical) gas parameter. The ex- 
perimental data points for /i(z/) are displayed in Fig. [51 
We observe a clear departure of the EoS from the mean- 
field prediction (dashed red line). At the largest mea- 
sured value of 1/ = 2.8 X 10~^ our data show a reduction of 
20 % of the pressure with respect to the mean-field result 
h^iy) = I'Kv'^ . In Fig.[2]we compare the measurements to 
the mean-field EoS, the LHY result (translated to the 
grand-canonical ensemble), and our QMC calculation. 
We observe that LHY theory accurately describes our 
experimental data and is hardly distinguishable from the 
QMC in the studied range of interaction strength, a point 
already put forward in a diffusion Monte Carlo simulation 
at even higher values of the gas parameter [2l| . We can 



3 



quantify the deviation of our data from mean- field theory 
by fitting the measured EoS with a function that includes 



a correction of order 



For this purpose we convert 



the energy E/N = 2TTh^ / {ma?)-na?' (l + a{na^Y/'^ 
to the grand-canonical EoS (see the supplementary ma- 
terial) and use q as a fit parameter in the resulting 
pressure P{^J^). The fit yields the value a = 4.5(7), 
which is in excellent agreement with the theoretical result 
128/(15V7F) « 4.81 in Eq. (P). Together with the mea- 
surement with composite bosons of this provides a 
striking check of the universality predicted by the expan- 
sion (jlj up to order vr 



22|. 



In the above interpretation we assumed that the zero- 
temperature regime has effectively been reached. To 
check this crucial assumption, we have performed finite- 
temperature path-integral Quantum Monte Carlo simu- 
lations 23| in the anisotropic harmonic trap geometry 
of the experiment with continuous space variables. The 
experimental atom number can be reached without dif- 
ficulty and pair interactions are described by a pseudo- 
potential. All the thermodynamic properties of the gas 
at finite temperature are obtained to high precision and 
without systematic errors. From a comparison of the ex- 
perimental density distributions with the QMC profiles 
at different temperatures (see FigH]), we find that the 
clouds in our experiments are colder than 0.25 Tc, where 
Tc is the condensation temperature of the ideal Bose gas. 
Thus thermal effects lead to an error in the EoS which is 
much smaller than the statistical error bars in Figl2l 

We now study the dynamics of a Bose gas subjected 
to a time-dependent interaction strength towards the 
strongly interacting regime. Besides its own interest, the 
investigation of this problem is important to assess the 
adiabaticity of the interaction sweep in the measurements 
described above. At strong interactions the enhanced 
three-body loss rate limits the practical duration of the 
sweep. In that violation of adiabaticity could lead 

to non-equilibrium density profiles that distort the mea- 
sured EoS. To address the question of adiabaticity, we in- 
vestigate the axial size of the Bose gas after a sweep of the 
scattering length. In Fig. [3] we plot the cloud size deter- 
mined by a Thomas-Fermi fit as a function of the sweep 
duration. The magnetic field is ramped approximately 
linearly in time, sweeping a/oo from an initial value of 
200 to different final values. Besides the experimental 
data we present theoretical results from a mean-field 
scaling solution 24, [2^ and from a solution of the hy- 
drodynamic equations incorporating the LHY EoS based 
on a variational scaling ansatz [26j . The latter shows 
a remarkable agreement with our experimental data for 
a < 3000 ao. For scattering lengths a/ao < 840 the radius 
is nearly constant for sweep durations TLUz/i^ir) > 1.5 
(r > 80 ms), indicating that the cloud follows the in- 
teraction strength adiabatically. For the largest value 
employed in the EoS study {a/ao — 2150), the atom loss 
is less than 4% and the cloud size after the r = 150 ms 



sweep {TUJz/{'2Tr) > 2.8) is 2.5 % smaller than the equilib- 
rium value. We have corrected for this systematic effect 
by rescaling the measured density uq for the determi- 
nation of the EoS, n = ri~^fio{riz) (with rj = 0.975 for 
a/ao = 2150). 
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FIG. 3: (Color online) Radius R of the Bose gas as a func- 
tion of the duration r of the interaction sweep. The radius 
R is normalized to the radius 7?* = aho(15A^A'^)^'^^ (where 
tsho ~ (Ji/mujzY^'^ and A = uJr/uJz)- N is the measured atom 
number at the end of each sweep. The final values of a/ao are 
380 (blue dots), 840 (purple squares), 2940 (red diamonds) 
and 4580 (green triangles). The solid (dashed) lines show the 
solution of a variational hydrodynamic approach (mean-field 
scaling solutions). The crosses show the predicted equilibrium 
beyond mean-field radii. 

The properties of the Bose gas for large values of na^ 
constitute a challenging open problem. Due to the exper- 
imental limitation imposed by three-body recombination, 
we access this region with a shorter sweep of duration 
tuJz/{2tt) = 1.35 (r ~ 75ms). In Fig. U] we plot the nor- 
malized radius of the Bose g function of the inverse 
scattering length aho/a. Deep in the mean-field regime 
(a < 800 ao) the ramp is adiabatic as the data matches 
the equilibrium Thomas-Fermi prediction. As the scat- 
tering length is increased, both non-adiabaticity and be- 
yond mean-field effects become important. A departure 
from the equilibrium result becomes evident above a scat- 
tering length of ~ 2000 ao. Taking into account the 
mean-field dynamics gives already an improved descrip- 
tion of our data (red dashed line). Even better agree- 
ment (up to values of a/ao ~ 5 000) is obtained with the 
variational approach incorporating the LHY correction 
as presented above (green solid line) p6| . 

Further insight is gained by considering data taken 
right on the Feshbach resonance and beyond. Because of 
the low densities of our samples, only half of the atoms 
are lost at the end of the sweep to the resonance (see 
squares in Fig. |4]). Universal thermodynamics at uni- 
tarity (where a diverges) has been conjectured for quan- 
tum gases [IB] and successfully checked experimentally 
for Fermi gases [slj. In the case of bosonic atoms the 
existence of a many-body universal state at unitarity is 
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FIG. 4: (Color online) Normalized cloud radius Rtf/R* 
(filled purple circles) and normalized atom number (open 
black squares) as a function of the inverse scattering length 
aho/a at the end of a 75 ms magnetic field sweep. The static 
mean-field prediction is plotted in solid red line, the mean- 
field scaling solution in dashed red, and the beyond mean- 
field scaling ansatz in solid green line. Inset: Zoom around 
the unitary limit. Predictions for the universal constant f are 
shown in up triangle 127], down triangle 28] and square [29l ]. 
The filled (empty) circles correspond to the radii normalized 
to the final (initial) atom number. The dotted black line is 
the linear interpolation at unitarity. 



still unknown. If universality proved to be correct, the 
only relevant length scale would be the interparticle spac- 
ing and the EoS would take the form fi cx ^n^/'^. 
Up to a numerical factor, this EoS is identical to that 
of an ideal Fermi gas and we can write ii = ^Ep (where 
Ef = h? /2m (67r^n)^/'^). When the scattering length is 
increased, the cloud is expected to grow in size. Due 
to the finite response time of the gas, it is reasonable 
to assume that the measured radius R is smaller than 
the equilibrium radius. From this inequality we deduce a 
lower bound for the value of ^ by interpolating our data 
at unitarity as presented in the inset of FigU) f > 0.44(8) 
33. This bound is satisfied for the predictions ^ — 0.66 
291 1 and for the upper bounds from variational calcula- 
and 2.93 



27]. 



tions, 0.80 [26 

Future work could focus on the measurement of the 
condensate fraction since the quantum depletion is ex- 
pected to be as large as ~ 8% for our most strongly 
interacting samples in equilibrium, and on finite temper- 
ature thermodynamic properties. Our measurements on 
resonance as well as future theoretical studies should give 
crucial insights on the unitary Bose gas. 
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[32] At unitarity, the cloud radius R would scale as A'^^''*'^^''^. 
The normalization R* used in Fig. 4 scales as N^^"' so 
that R/R* (X ^i/4jY-i/3() 

very slowly depends on atom 
number. In order to take into account the changing atom 
number near unitarity and obtain a safe experimental 



lower bound on ^ oc {R/ R*)'^N'^^^^ , we minimize both 
R/R* and N'^'^^\ This is done by taking for R/R* the 
initial atom number (empty circles in the inset of Fig.4), 
and the final, for A^^/is 



